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Abstract. The LM method (Luri, Mennessier et al 
|1996| ), designed to exploit the Hipparcos data to obtain 
luminosity calibrations, is applied to derive luminosity 
calibrations for RR Lyrae and classical Cepheids. From 



Lyrae 18. 31™ (Fernley et al., 1998) (direct determination), 
18.63™ dGratton et al., 1997D, 18 65"' (|Reid, 1997]) (indi- 



rect determinations obtained from subdwarf-sequence fit- 
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these calibrations the distance to the Large Magellanic 18.44 — 18.57™ ( 


Madore & Freedman, 1998 


), 18.72™ 


(Pa- 


Cloud 'LMC) is estimated. The distance moduli provided turel et al., 1997 


), 18.70™ ( 


Feast & Catchpole, 1997 


)■ 



by the two calibrations are in good agreement, giving a 
value of ~ 18.3™, while several previous calibrations us- 
ing Hipparcos data provided inconsistent results between 
both types of stars. This result suggest that the Hubble 
constant should have a value of Hq ~ 79 km s~^Mpc~^. 
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In this paper luminosity calibrations for both RR 
Lyrae and classical Cepheids are obtained using the LM 
method applied to Hipparcos data. The results provide 
compatible values for the LMC distance modulus. 



2. The LM method 

The LM method (Luri, Mennessier et al., 1996) is based 



1. Introduction 

The calibration of the absolute magnitudes of RR Lyrae 
and the classical Cepheids is the first step in the deter- 
mination of the extragalactic distance scale, and the re- 
cently released Hipparcos data ( ESA, 1997 ) allow, for the 
first time, its determination on the basis of trigonomet- 
ric parallaxes. However, in spite of the high accuracy of 
these data, few of these stars have precise trigonometric 
parallax measurements: only 12 RR Lyrae and 6 classi- 
cal Cepheids have relative errors in trigonometric paral- 
lax smaller than 30%. Due to this limitation, and to the 
intrinsic difficulty of determining distances and absolute 
magnitudes from trigonometric parallaxes (several biases 
may arise from the effects of the observational errors and 



on the Maximum-Likelihood estimation. It includes a de- 
tailed model of the luminosity, kinematics and spatial dis- 
tribution of the sample and takes into account its observa- 
tional censorship and observational errors, thus providing 



sample censorship, see Brown et al. (1997)), a careful sta- 
tistical treatment of the data is required to obtain reliable 
calibrations. 

The difficulty of these estimations is illustrated by 
the wide range of values for the distance modulus of the 
Large Magellanic Cloud (LMC) obtained from published 
luminosity calibrations using Hipparcos data: from RR 



estimations free of biases due to these two factors (Luri 
|fc Arenou, 1997 ). The interstellar absorption is taken into 
account by using the Arenou et al. (1992) 3D model. 

Using the LM method, the parameters of the model 
used are estimated. The estimation uses all the available 
information for the stars in the sample: apparent magni- 
tude, galactic coordinates, trigonometric parallax, proper 
motions, radial velocity and any other relevant parameter 
such as metallicity or period. The use of all the observa- 
tional data is specially important in the present case be- 
cause parallaxes alone would not provide a precise enough 
calibration (their relative errors being high, even with the 
Hipparcos high-precision astrometry). Furthermore, as the 
estimation is done by Maximum-Likelihood, the informa- 
tion given by these observational data is included through 
the Probability Density Function (PDF) defined by the 
model and the observational errors. Consequently, each 
individual piece of data has its own "intrinsic weight" in 
the solution ^ and there is no need, as in other methods 
used for absolute magnitude calibration, for any external 
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^ the relative contribution of parallaxes and proper motions 
to our solutions will be assessed in future papers 
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system to weight the contribution of, say, parallaxes or 
proper motions in the estimation. 



3. RR Lyrae 



of the Halo RR-Lyrae luminosities from two different 
methods. After averaging them they adopt a value of 
< My >= 0.77 ±0.15 at [Fe/H] = -1.53 is adopted. The 
differences with our results can be accounted for by the 
different criteria used to separate Halo and Disk. While 
Fernley et al. (199S) use an a priori metallicity criterion 
to divide the sample into Halo and Disk, our separation 



The LM method was adapted to determine a mean ab- 
solute magnitude and the corresponding dispersion cta/ 
for the RR Lyrae stars. The distribution of metallicities 
was modeled (and fitted) using normal distributions. To 
model the kinematics of the sample a velocity ellipsoid 
with means (Uq, Vb, Wo) and dispersions (cr;/, cry, ctw) was 

adnptP d An PvpmiPntia.1 gfllartir Hkk with sprIp hpight Zp ties using Hipparcos data a re inconsistent with ours ( |Rcid 



is part of the fit, taking simultaneously into account the 
luminosity, the kinematics and the metallicity of the stars. 

Other recent estimates for the halo RR Lyrae luminosi- 



was usl^d to describe the spatial distribution 



On the other hand, the apparent magnitude selection 
of the sample was also taken into account. The Hipparcos 
catalogue was designed to be complete up to an apparent 
magnitude varying on galactic latitude and spectral 
type, and for fainter magnitudes very heterogeneous 
selection criteria were used. In the case of RR-Lyrae 
the criteria used to complete the catalogue up to the 
Hipparc os m agnitude limit is described in Mennessier & 
Baglin ( 1988| ) and, furthermore, six previously unknown 
RR-Lyrae were found. This observational censorship was 
modeled in the LM method by assuming the sample to be 
complete up to an apparent magnitude Vc (determined 
at the same time than the rest of the parameters) and 
with a linear decrease in completeness up to the apparent 
magnitude limit, reflecting the fact that fainter RR-Lyrae 
have a smaller probability to be included. 

The data used for the RR Lyrae calibration comes from 
two so urces: astro metric data from the Hipparcos Cata- 
logue ( ESA, 199^ ) and intensity- mean V apparent magni- 
tudes (calculated from the Hipparcos data), metallicities 
and radial velocities from the compilation of Fernley et 
al. ( 1998 ). There are 186 RR Lyrae stars in the Hipparcos 
catalogue, 6 of them newly discovered. The Fernley compi- 
lation contains 144 stars (125 RRab and 19 RRc) reliably 
classified as RR-Lyrae, which constitute our sample. 

The LM method identified two main groups, consti- 
tuting the 91% of the sample. The first group corresponds 
to the Halo population and the second to the Disk pop- 
ulation. The mean magnitudes and metallicities for these 
groups are listed in Table 0. 

Table 1. Mean absolute magnitudes and metallicities for 
RR-Lyrae 



< Mv > 



[Fe/H] 



% of the sample 



Halo 
Disk 



0.65 ±0.23 
0.13 ±0.49 



-1.51 ±0.06 
-0.45 ±0.07 



78.3 ± 2.4 
12.7 ± 1.6 



Our results can be c ompared with those reported 
by Fernley et al. (|l998| ). They obtain an estimation 



1997; Gratton et al., 1997), giving brighter mean absolute 



magnitudes. However, they are indirect estimates based 
on determinations of the subdwarf sequence and they in- 
clude a posteriori corrections of parallax biases that can 
degrade their precision ( Brown ct al., 1997 ). 



4. Classical Cepheids 

For these stars we consider a period-luminosity (PL) rela- 
tion (Eq. 0): 



<AF>=A + B log(P) 



(1) 



It was assumed that for each value of the period the 
individual values of < My > are distributed normally 
around the PL relation with a dispersion aM ■ The periods 
were modeled (and fitted) using normal distributions. The 
kinematics, spatial distribution and apparent magnitude 
selection were modeled as explained in Section ^. The val- 
ues of Oort's constants and the Sun's galactocentric dis- 
tance were not determined but adopted to be A = 14.4, 
B = —12.8 km s~^kpc~^ and Rq = 8.5 kpc. 

The sample was formed by selecting the classical 
Ceph eids ((5-Cepheids) of the Hipparcos catalogue ( ESAj 
1997 ). The known sinusoidal (5-Cepheids (overtone 
Cepheids) were eliminated. All data (including periods) 
were taken from the Hipparcos catalogue except the radial 
velocities, taken from the Hipparcos Input Catalogue 
(Turon et al., 1992). The arithmetic-mean apparent 



magnitudes given by Hipparcos were compared with 
the intensity-mean apparent magnitudes given in the 
David Dunlap Observatory Database of Galactic Classical 
Cepheids ^ and no systematic difference was found (mean 
difference 0.01™ ± 0.01). Thus, the Hipparcos data were 
preferred due to their higher homogeneity. The final 
sample contains 219 stars. 

Two determinations of the PL relation were obtained: 

Relation 1: following the approach taken by Feast & 
Catchpolc (1997), the PL slope (B) was fixed to the 



value for the LMC, B = -2.81 (Caldwell fc Laney 



1991). The underlying hypothesis is that the slope of 
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the PL relation is (except for a small metallicity cor- 
rection) universal, so the slope for the LMC Cepheids 
can be used and only the zero point of the relation 
remains to be determined. 
Relation 2: both the slope and the zero point are deter- 
mined. 

In both cases the LM method identified a small secondary 
group, but the most part of the sample (91%) belongs to 
the main group. The two solutions obtained for the PL 
relation of this main group are presented in Table ||. 

Table 2. Period-luminosity relations for the classical 
Cepheids 

Relation 1 



< Mv >= -2.81 log(P) - (1.05 ± 0.17) 

Relation 2 

< Mv >= -2.121og(P) - 1.73 
e<M„> = 0.20 0.08 log(P) 



In the case of Relation 2, the slope and zero point of 
the PL relation were not used as parameters directly deter- 
mined by the method due to the high correlation between 
them could degrade the precision of the numerical method 
used to maximize the likelihood. Instead, two points of the 
PL relation (at two arbitrary values of the period) were 
determined (thus defining the linear relationship) and the 
slope and zero point were calculated from them. Conse- 
quently, the errors in the estimates of the zero point and 
the slope cannot be given independently and, instead, an 
estimation of the expected error in the absolute magnitude 
(e<M„>) is given as a function of log(P). 

For the Cepheids, unlike the RR Lyrae, the err ors in 
the interstellar absorption from the Arenou et al. ( 1992 ) 
model (hereafter AGG) can be high (most of the stars are 
located in the galactic plane and at higher distances than 
the RR Lyrae). To obtain the value of the intrinsic dis- 
persion we should take into account that the value of the 
magnitude dispersion given by the LM method is the re- 
sult of this dispersion aM and the errors in the estimation 
of the interstellar absorption ■ <^m total — '^m + '''a — 
(0.8™ ±0.1)2. 

The AGG model provides estimations of the errors in 
the values of the interstellar absorption. Using these es- 
timations to correct the total dispersion, the value of the 
dispersion of the sample around the PL relation can be es- 
timated as (Jm = 0.4™ ± 0.2. In any case, the PL relations 
obtained do not depend on the value of this parameter, as 
shown by Monte-Carlo simulations. 

A recent result for the PL rela tion fr om Hipparcos data 
is the one of Feast & Catchpole ( 1997 ) (hereafter FC): 



< My >^ -2.81 log(P) + (-1.43 ± 0.10) 

This result can be compared with our Relation 1 (both 
rely on the hypothesis of a known slope B = —2.81). Our 
zero point is 0.38™ fainter than that given by FC but, 
before a discussion of this difference some details about the 
FC approach are necessary. To determine the zero point of 
the PL relation FC use the following method. Given Eq. 
and Po gson's law, the following relation holds: 



10 



0.2A 



0.0l7rl0°-2[<v'>o-siog(p)] 



(2) 



where < V >o is the intrinsic apparent magnitude, i.e. cor- 
rected for interstellar absorption. For each star the quan- 
tity Q = lO'^'^'^ can be estimated and the zero point of the 
PL relation A calculated from the mean value obtained for 
all the stars. 

This method of estimating the zero point of the PL re- 
lation has the advantage of avoiding the direct calculation 
of absolute magnitudes from parallaxes, which can lead to 
a bias (even when using Hipparcos unbiased parallaxes) 
if not treated properly (Brown et al., 1997). Instead, the 
parallaxes are directly averaged and the zero point esti- 
mated from the average, minimizing this source of bias. 
However, the method is highly sensitive to any error in 
the exponent of the right hand side of Eq. |^, including 
any effects on the magnitude distribution (like Malmquist 
bias) or the reddening correction. On the other hand, the 
weighting system used by FC to obtain the mean value of 
Q for the sample can have some undesired side-effects: as 
the weigth of each star is proportional to ^ , being Qi the 
individual value of Q for the star, stars with low walues 
of Qi are favoured in the final solution. Furthermore, due 
to the weighting only a (relatively small) fraction of the 
sample significantly contributes to the solution, so arising 
the issue of how representative of the whole population 
are these contributing stars. 

The impact of these effects on FC method is difficult 
to evaluate, but Monte-Carlo simulations of realistic 
samples show that the zero point given by FC could be 
slightly (about 0.05 — 0.1™) too bright due to them, con- 
tributing to explain in part the difference with our results. 

On the other hand, when the LM method is applied 
using the absorption correction method given in FC 
instead of using the Arenou et al. ( |1992D model, a 
value of a ^ total 

= 0.7™ ± 0.2 is obtained. This result 
suggests that the combination of gm and errors in the 
absorption estimation ua^ fc gives a total dispersion 
higher than the estimated by FC. The PL relation, 
< Mv >= -2.041og(P) - 1.74, and the kinematics and 
scale height obtained do not differ significantly from the 
results in Relation 2. 

Our second relation (Table ||) gives a slope of the PL 
relation less steep than the one given by Caldwell & Laney 
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( |l99l| ), but consistent with the results of Szabados (1997) 
from Hipparcos data for nine non-binary Cepheids with 
short periods. 

Further analysis to determine the slope of the PL re- 
lation are being carried on using the LM method and 
the preliminary results suggest a different behavior in the 
short and long period regions, possibly due to the effects of 
undetected overtone cepheids in the short period region. 

5. The LMC distance modulus 



The calibrations presented in this paper were used to de- 
termine the mean distance modulus of the LMC. The re- 
sults are presented in Table ^ and they were obtained as 
follows: 

RR Lyrae: to calculate the distance modulus of the LMC 
using RR-Lyrae data, a value of the slope of the 
metallicity- luminosity relation is needed. Although 
the value of this slope could be determined using the 
LM method, an adopted value was used here, leaving 
for a forthcoming longer paper the discussion of this 
parameter. Notice, however, that the mean magnitude 
determined here corresponds to a value of metallicity 
(—1.51) close to the mean value of the LMC RR-Lyrae 
(—1.8) so the resulting distance modulus does not de- 
pend strongly on the value of the slope adopted. 
Following the approach of Fcrnley et al. (|1998| ) a slope 
of 0.18 was adopted. Using this value and the results 
for the Halo RR-Lyrae given in Table |l|, a metallicity- 
luminosity relation was obtained and applied to the 
RR Lyrae data given in Walker (1992) (individual red- 
dening estimates used). 

Cepheids (FC revised): the FC estimation of the LMC 
distance modulus was changed by 0.38™ to reflect the 
change in zero point in our Relation 1. 

Cepheids (Rel. 1 & 2): the PL relations given in Table I 
were applied to the Cepheid data given in Paturel et 
al. (|1997 ); a mean reddening correction of Eb-v = 0.1 
(Freedman et al., 1994) and a metallicity correction of 



-1-0.042 (Laney fc Strobie, 1994) were applied 



Table 3. Distance modulus of the LMC using this paper's 
luminosity calibrations 



M 



RR Lyrae 

Cepheids (FC revised) 
Cepheids (Rel. 1) 
Cepheids (Rel. 2) 



18.37 ±0.23 
18.32 ±0.17 
18.29 ±0.17 
18.21 ±0.20 



The results of this paper reconcile the distance modu- 
lus estimations of the LMC based on RR Lyrae and those 
based on the classical Cepheids. Moreover, they are consis- 
tent with the upper limit of 18.44 ±0.05 derived by Gould 



k Uza ( |1997D from the analysis of the SN 1987A super- 
nova "light echo" . The adoption of a value of 18.3™ for the 
distance modulus implies that the Hubble constant should 
now have a value of = 79 km s~^Mpc~^, in contrast 
to the value of iJo — 73 km s~^ given by Freedman et al. 



(1997) 
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